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Dimensional recurrence relations: 

an easy way to evaluate higher orders of expansion in e 
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Applications of a method recently suggested by one of the authors (R.L.) are presented. This method is based 
on the use of dimensional recurrence relations and analytic properties of Feynman integrals as functions of the 
parameter of dimensional regularization, d. The method was used to obtain analytical expressions for two missing 
constants in the e-expansion of the most complicated master integrals contributing to the three-loop massless 
quark and gluon form factors and thereby present the form factors in a completely analytic form. To illustrate its 
power we present, at transcendentality weight seven, the next order of the e-expansion of one of the corresponding 
most complicated master integrals. As a further application, we present three previously unknown terms of the 
expansion in e of the three-loop non-planar massless propagator diagram. Only multiple ( values at integer points 
are present in our result. 



1. Introduction 

Quite recently a new method pQ of evaluating 
Feynman integrals was suggested. It is based on 
the use of dimensional recurrence relations [2] and 
analytic properties of Feynman integrals as func- 
tions of the parameter of dimensional regulariza- 
tion, d. Here we are going to describe new results 
of its application. In the next section we give a 
summary of the evaluation of the two previously 
analytically unknown coefficients of transcenden- 
tality weight six in the e-expansion of two master 
integrals contributing to the three-loop massless 
form factors [3]. We also present, at transcen- 
dentality weight seven, the next order of the e- 
expansion of one of the corresponding most com- 
plicated master integrals. In Section 3, we report 
on the evaluation of higher orders of expansion 
in e of the three-loop non-planar massless prop- 
agator diagram and, in Section 4, we discuss our 
results. 
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2. Evaluating master integrals for 
three-loop form factors 

The integration-by-part reduction reduces the 
problem to the calculation of a small number 
of master integrals. All the master integrals 
apart from three most complicated master inte- 
grals contributing to the three-loop massless form 
factors have been analytically evaluated in [415] . 
About one year ago, one of the three most com- 
plicated master integrals (called A 9j i in |4l5|6|7j ) 
and the pole parts of Ag^ and Ag^ shown in 
Figs. 1 and 2 were evaluated analytically, while 
the e° parts of Ag^ and Ag y2 were evaluated nu- 
merically — see |6l7j . 

In [3], the two missing ingredients, i.e. the fi- 
nite parts of Ag t 4 and Ag t2 , were evaluated ana- 
lytically. According to the method of [1] it is nec- 
essary, before evaluating Ag^ and Ag t2 , to know 
all lower master integrals. They are shown in the 
same figures. Four rows of diagrams in each fig- 
ure correspond to complexity levels 0, 1, 2 and 

3. Details of the calculation can be found in pQ. 





Figure 2. Master integrals for A 
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Figure 1. Master integrals for Ag A . 



Here are the corresponding results: 
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An independent calculation of the form factors 
was performed quite recently [5] and the agree- 
ment with the previous results was established. 
Motivated by a future four-loop calculation the 
authors calculated also the subleading 0(e) terms 
for the fermion-loop type contributions. We think 
that the analytic calculation of the whole 0(e) 
part of the form factors is feasible at the moment. 
To illustrate this possibility we have calculated 
the 0(e 2 ) order of the £-expansion of one of the 



corresponding most complicated master integrals, 
Agi, which is of transcendentality weight seven: 



A 9A (4-2e) = e 
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3. Non-planar massless propagator dia- 
gram 

It is also natural to try to apply the same 
method to evaluate the three previously analyt- 
ically unknown coefficients of transcendentality 
weight six in the £-expansion of three master inte- 
grals contributing to the three-loop static quark 
potential |9I10I11I12I13TT4] . These are Jii,Ii 6 , ha 
in Fig. 3. 

This work is in progress [15] . Here we would 
like to present, as a by-product of this activity, 
new results for one of the master integrals which 
are lower than Iig. Let us consider In which is a 
master integral for three-loop massless propaga- 




Figure 3. Most complicated master integrals for 
the three-loop static quark potential. 



tor integrals. The famous value 20£(5) at d — 4 is 
known for a long time |16j . The linear term in the 
e-expansion was obtained in [17]. The quadratic 
term was recently evaluated [18 . To illustrate the 
power of the present method we have evaluated 
the three next terms so that we have a result up 
toe 5 . 

The diagram itself and the corresponding four 
lower diagrams are shown in Fig. [4] The labelling 






Figure 4. Master integrals for 7 17 . 



for Iii = Ng 3 and its lower master integrals is 
taken from the future publication [15] . This mas- 
ter integral itself has complexity level 2, the mas- 
ter integral P§^ has complexity level 1, while the 
other three master integrals can be expressed in 
terms of gamma functions at general d. 

The lowering recurrence relation for this inte- 
gral provides the following expression of AQ 3 
in terms of master integrals in d dimensions: 
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We obtain the following result to order e 5 : 
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where Q(n,rn) are multiple zeta values (see, e.g., 
19 ). Observe that the pulling out the standard 
prefactor 1/(1 — 2e) in results for massless propa- 
gator integrals (see, e.g., [16I17I20J ) provides the 
homogeneous transcendcntality weight in all the 
orders of the expansion in e. We see this property 
in our new piece of the result, i.e. in the e l terms 
with i — 3, 4, 5. Let us emphasize that this prop- 
erty is very useful when using PSLQ [21] because 
the number of constants that can be present in a 
result is essentially reduced. In fact, we indeed 
arrived at the above result using PSLQ but could 
here proceed even without the homogeneous tran- 
scendentality weight because, within the method 
of [T], we obtained, for J17, a very well convergent 
double series so that we could obtain the accuracy 
of thousand of digits or more. For the same rea- 
son, we can obtain higher terms of the expansion 
in e, i.e., e 6 etc. 

If one pulls out a more sophisticated prefactor 
[16120] (see the next equation) the pure 7r 2 factors 
will not be present in the result. We see that this 
property is satisfied indeed: 
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Let us observe that in [22] it was proven that 
the coefficients in the ^-expansion of planar mass- 
less propagator diagrams up to five loops should 
be expressed in terms of multiple zeta values, 
while the non-planar graphs may contain, in addi- 
tion, multiple sums with 6th roots of unity. How- 



ever, in the above result, only multiple zeta values 
are present. 

4. Conclusion 

Let us emphasize that in the present paper 
we used method of Ref.pQ in combination with 
several other methods. We already mentioned 
PSLQ. Then, as it was explained in |1I3] . it is im- 
plied that the IBP reduction is solved for a given 
problem so that one can obtain dimensional re- 
currence relations. To do this we used the code 
called FIRE gl] and the code based on ^5\ but, 
of course, one can use other methods. Moreover, 
we use a sector decomposition [26 27 28^ imple- 
mented in the code FIESTA |28)29j to determine 
the position and the order of the poles in the basic 
stripe. To fix remaining constants in the homoge- 
nous solution of dimension recurrence relations 
we apply the method of Mellin-Barnes represen- 
tation j3QI31|32j . 

The method used in the present calculations 
looks quite promising and we expect it to be ap- 
plied not only in the problems discussed above 
but also in many other situations. 
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